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Abstract. The aim of electrical impedance tomography is to form an image 
of the conductivity distribution inside an unknown body using electric bound- 
ary measurements. The computation of the image from measurement data is 
a non-linear ill-posed inverse problem and calls for a special regularized algo- 
rithm. One such algorithm, the so-called D-bar method, is improved in this 
work by introducing new computational steps that remove the so far necessary 
requirement that the conductivity should be constant near the boundary. The 
numerical experiments presented suggest two conclusions. First, for most con- 
ductivities arising in medical imaging, it seems the previous approach of using 
a best possible constant near the boundary is sufficient. Second, for conduc- 
tivities that have high contrast features at the boundary, the new approach 
produces reconstructions with smaller quantitative error and with better visual 
quality. 



1. Introduction 

The aim of electrical impedance tomography (eit) is to form an image of the conduc- 
tivity distribution inside an unknown body using electric boundary measurements. 
Applications of eit include medical imaging, nondestructive testing and subsurface 
monitoring. See [T6^ for an overview of eit. The computation of the image from 
measurement data is a non-linear ill-posed inverse problem and calls for a special 
regularized algorithm. In this work we improve one such algorithm, the so-called 
D-bar method, by removing the so far necessary requirement that the conductivity 
should be constant near the boundary. 

The mathematical model behind eit is the inverse conductivity problem intro- 
duced by Calderon in pj: . We discuss here the inverse conductivity problem in the 
following two-dimensional form: let Qi = 1^(0, ri) C be the disc with center at 
origin and radius ri > and consider a strictly positive, real- valued conductivity 
a e C^(r^i). Maintaining a voltage distribution / on the boundary dfti creates a 
voltage potential u that solves the following Dirichlet problem: 

. . J V • {aVu) = in Qi, 

\ u = f on dQi. 

The resulting distribution of current through the boundary is 
(2) Aaf = cT—\dn, 

where u is the outward unit normal and A^r is the Dirichlet-to-Neumann (dn) map. 
Calderon's problem is to reconstruct a from the knowledge of Act. 

Many numerical methods have been suggested in the literature for the recon- 
struction of a in the above setting. In this work we concentrate on the so-called 
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D-bar method. Other approaches include hnearization [5j [H [17], iterative regu- 
larization [20l EH |37l [26] , statistical inversion [25l |45] , resistor network methods 
pT| . convexification [4T, layer stripping |49l |50] and Teichmiiller space methods 
|3ll [35l [36] . Also, there is a large body of work concentrating on recovering partial 
information on <j, see [9] |10] for a survey. 

Theoretical foundation of the D-bar method for eit in dimension two was intro- 
duced by Nachman in [42 , where a constructive proof for recovering cr G W'^'^iyti) 
from Act was given for p > 1. The result was later sharpened by Brown and 
Uhlmann [13 to cover W^'^{Vti) conductivities with q > 2; the proof was aug- 
mented with constructive steps by Knudsen and Tamasan in [29]. Finally, Astala 
and Paivarinta answered Calderon's question in its original form by describing a 
constructive procedure for recovering a G L'^i^Qi) in [3, 4 ; numerical implementa- 
tion of this approach is described in ^ 2^. Thus there are several variants of the 
D-bar method for two-dimensional EIT. In dimension three, the theory of D-bar 
reconstruction is developed in [191 [SI HH Ell E] • 

The above theoretical results on the D-bar method assume the knowledge of the 
infinite-precision data Act- However, the starting point of practical inversion is a 
finite-dimensional and noisy approximation L^. to A^r. Since the eit problem is 
severely ill-posed, or sensitive to measurement noise, any practical reconstruction 
method needs to be robust against errors in measurement data. The first robust 
D-bar algorithm (based on [42 ) was given in [48 , and it has been refined and 
analysed in [38, ,301 [SH [32 • The method has been successfully tested on a chest 
phantom in [23] and on in vivo human chest data in [24] . Numerical D-bar method 
based on [T3", "29] was reported in [28 . The above methods are two-dimensional; 
three-dimensional computations are described in [6, 12 . Robustness is ensured in 
all of these practical D-bar reconstruction methods by truncating scattering data, 
a step that can be viewed as nonlinear low-pass filtering. 

In addition to being robust against noise, a reliable eit algorithm needs a reg- 
ularization analysis. Such an analysis is provided for the two-dimensional D-bar 
method in [33j, where an explicit formula is given for choosing the truncation ra- 
dius as function of noise level. 

Practical D-bar methods have been until now implemented by fitting an optimal 
constant to the possibly nonconstant trace cr|afii. Quite good results have been 
obtained both with laboratory data [23j and in vivo patient data [24J. However, 
in applications exhibiting large conductivity changes near or at the boundary the 
constant-fitting approach may not be good enough. Our aim here is to remove the 
assumption "a = 1 near the boundary" from the two-dimensional D-bar algorithm 
based on [42j using an additional procedure we call boundary correction. 

Let us review the infinite-precision boundary correction procedure given in [42] . 
The starting point is the DN map A^r of a conductivity a G Take r2 > ri 

and set Q2 = ^(0,r2). The conductivity a is extended outside Qi by 



where we can choose any a G W^'^(l]2 \ ^1) with the properties = cr|aOi and 

{da / diy)\dn^ = {da / diy)\dni and a = 1 near 9(^2- This way 7 G W'^'^{ft2) whenever 



(3) 




cr(x), when x G 

when x G ^^2 \ ^1, 
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(J G W'^'^iyti). Define two Dirichlet problems: 



The boundary corrected D-bar method for cr G W'^'^(yti)^ assuming infinite-precision 
data, is described as fohows in [42, Section 6]: 

(a) Reconstruction at the boundary. Recover the trace cr\d^-^ and 
the normal derivative {dcj / dv)\dn^ from A^r- 

(b) Extension of conductivity. Using (a) and ([3]), extend the conductivity 
to 7 G W'^'^{Vt2) such that 7 = 1 near 80.2 and m.ixen2 ^ c > 0; 

(c) Calculation of outer DN map. Write the DN map of 7 G W^'^(l]2) 
in terms of A^r using (j6j); 

(d) Reconstruction using the D-bar method. Reconstruct 7 G W'^'^{VL2) 

from the infinite-precision data A^ following [42^ . 

The practical starting point of reconstruction is the noisy approximate data L^, 
and the Steps (a-d) above cannot be directly followed. We suggest the following 
robust procedure for boundary correction: 

(a^) Approximate reconstruction at the boundary. Recover numerically 
a smooth function g G (7^(91^1 ) with the property g ~ cf\q^^ as explained 
in [43]. Omit recovering {dcr / dv)\dni since it is an unstable step [43]; 



(b^) Simple extension of conductivity. Construct a G C^(1^2 \ ^1) 
satisfying m^^^^^^^^{x) > c > and cf\q^^ = g and a = 1 near 91^2- 
Use a in ([3]) to extend the conductivity to 7 G L^(r^2); 

(c^) Approximate calculation of outer DN map. Write approximate 
DN map in terms of L^r using a matrix approximation to (j6|); 

(d^) Reconstruction using regularized D-bar method. Reconstruct 7 
from using the regularized D-bar method described in [33J. 

The main concern about the procedure (a'-d^) is that the extension of <j to 7 will be 
in general discontinuous at dOi, and thus 7 violates the assumptions of the D-bar 
method used in (d^). However, there is both theoretical and experimental evidence 
suggesting that the step (d^) should give reasonable results even in this case [23l 




(5) A'Vi='^^U, i,i = l,2. 
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[24t [32| [3T] . Another potential problem arises from the inverse operator in formula 
dni, as the proof of invertibility [42l Proposition 6.1] in the extended conductivity 
produced by step (b^). One possibility would be to use [27, Lemma 2.1.3.] instead 
of (|6j) as the basis of step (b'). However, in our computational experiments the use 
of (jni seems not to be a problem even in the case of discontinuous conductivity 
extensions. 

This paper should be viewed as a report of computational experiments suggesting 
the practical usefulness of the boundary correction step in applications where the 
conductivity varies strongly near the boundary. Hopefully the computational results 
presented below will act as motivation for further theoretical study of practical 
imaging algorithms for EIT. 

We remark that the boundary correction method is applied in this paper only 
in the case of l^i being a disc. This is not a serious lack of generality, though: we 
presume that other domains than discs could be treated combining the methods 
described in [39", '40^ with the boundary correction. 

This paper is organized as follows. We present our method of simulating contin- 
uum model EIT data in Section O The details of implementation of Steps (a') and 
(b^) and (c^) are discussed in Sections [3] and H] and [5l respectively. A brief outline 
of the regularized D-bar method is given in Section [6l Our practical boundary cor- 
rection method is illustrated by numerical examples in Section [71 and we conclude 
our results in Section [51 



2. Simulation of measurement data 

Let lZ(j : H~^^'^{dVti) H^^'^{dVti) denote the Neumann-to-Dirichlet map of a, 
where spaces consist of functions with mean value zero. We have IZaQ = 
u\dni^ where u is the unique H^{Qi) solution of the Neumann problem 

Ou 

V • a\/u = in 1^1, 7—— = ^ on dQi^ 

ov 

satisfying J^^^ uds = 0. We note two key equalities concerning A^r and T^cr- Define 
a projection operator P(j) := (j). Then for any / G H^^'^{dVti) we have 

PKf = \dni\-^ Jg^^a^ = V • ^^ii = 0, so actually A^ : H^^^{d^i) 
H~^^'^{dfti). From the definitions of A^r and 71a- we now have 

(7) AaUa = I : H-^^\dni) ^ H-^^^dQi), 

(8) n^Aa = I-P :H^^\dni)^H^^\dni). 

Given a and A/" > 0, we define a matrix : C^^ as follows. We use 

a truncated orthonormal trigonometric basis for representing functions defined at 
the boundary dQj: 

(9) ^(")(^) = ^^e»^ n = -N,...,N, j = l,2. 

Note that J^^ 4>^j^^ds = for n 7^ 0. Then solve the Neumann problem 

(10) V-ctVw^"^ =0 in (T^^ = </.$"^ on aS^i, 

ou 
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with the constraint Jq^^ u(^^ ds = 0. Define Ra^ = [u{£^n)] by 

(11) u{i,n)= [ u^l'Ui^ds. 

Here i is the row index and n is the column index. 

The matrix represents the operator IZa approximately. We add simulated 
measurement noise by defining 

(12) R', := Ra + cE, 

where £^ is a 2N x 2N matrix with random entries independently distributed ac- 
cording to the Gaussian normal density A/'(0, 1). The constant c > is adjusted so 
that \\Rl — Ri\\/\\Ri\\, where ||-|| is the standard matrix norm and Ri is the ND 
-map for the unit conductivity, is greater than the relative error caused by FEM and 
of the same order of magnitude as 0.0017% (signal to noise -ratio of 95.5 dB) , the 
noise level of the ACTS impedance tomography imager of Rensselaer Polytechnic 
Institute [18j. 

We can now easily compute the corresponding noisy matrix representation L^. 
for the DN map Act. Namely, define 

then L J is a matrix of size 2 A x 2 A. We should add appropriate mapping properties 
for constant basis functions at the boundary according to the facts 



A,1 = 0, / A,fds = 0. 

This is achieved simply by setting (in Matlab notation) 

L^(1:A, 1:A) £i(l:A, (A + l):end) 



(13) K'r- 







Lj((A + l):end,l:A) Lj((A + l):end, (A + l):end) . 
where the zero block matrices above have various (but obvious) sizes. 

3. Approximate reconstruction at the boundary 

The trace a\QQ^ can be approximately reconstructed in the following way 
Define hM,(3{0) = e'^^r]{0 - /3), where 



(14) rj{0) = 



d{tzO - 7r/2)"(A^(9 + 7r/2)" cos(/^l9), for - 7T/{2tz) <0 < 7i/{2k), 
0, otherwise 



is a non- negative cut-off function satisfying j^^^ rf{0)dO = 1. Now the mollified 
trace (cr7^^)|aJ7i (Z^) can be calculated with 

(15) / ar]'^ds = lim / hM,(3KhM,(3ds. 

We get the approximation g ^ cr\QQ^ by calculating ([15]) with different angles P and 
using a finite M in the right side of ([T5]) . 
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Another approach to reconstructing (j\dni is the layer stripping method intro- 
duced in [49]. 

4. Simple extension of conductivity 

The starting point here is a given approximation g : dVti ^ R to the trace cf\q^^ of 
the conductivity a G (7^(1^1 ) of interest. The aim is to construct a strictly positive 
conductivity a : 1^2 \ satisfying d\drLi — 9 and a = 1 near the outer boundary 
91^2, and then use formula ([3]) to define 7. 

We extend a to 7 using the following extension in polar coordinates: 



(16) 7(P,^) 



a(p,6>), P<ri, 
{g{0)-l)fm{p)^h n<p<re, 

^1, Te < p< r2, 

where ri < rg < r2 is some radius and fm{p) > is a suitable third-degree poly- 
nomial satisfying /m(^i) = 1 and fm{re) = 0. Note that 7 is twice continuously 
differentiable apart from possible discontinuity at p = ri , and equals constant 1 in 
the annulus < p < r2- 

5. Approximate calculation of outer dn map 

Using the functions (j9j), a given function / : dVti C can be approximately 
represented by the vector 

/ = m-N), n-N + 1), . . . , m - 1), f{N)f, /(n) = / f<p^Us, 



and the DN maps A*-^ can be approximated by the matrices L*-^ = [gij{^^^)] with 



(17) g^J{^.n)= I a^lonA'Us, 



where u^^-' denotes the solution to dH with u^^-'Iqq. = (t>f^^ - Again I is the row 
index and n is the column index. Now the matrix can be calculated by 

(18) Ll^ = L^^ + L^\L%-L^')-^L^\ 

provided that the matrix — L^^ is invertible. Formula ([18]) is a finite-dimensional 
approximation to (|6j) . 

6. Regularized D-bar method 

In this section we explain how to reconstruct a conductivity 7 in a regularized way 
from a noisy measurement matrix under the assumptions 7 G C^(1^2) and 7 = 1 
in a neighborhood of 80.2- 

If we had the infinite-precision data A^ at our disposal, we could follow the re- 
construction procedure in [42 . First we would solve the boundary integral equation 

(19) ^{ • , k)\Qn, = e''^ - Sk{A^ - A,)^{ • , k)\on,^ 

in the Sobolev space H^^'^{dVt2) for ah G C \ {0}. In formula ([IHl), 5'/^ is a 
single-layer operator 

{Sk(t)){x) \= / Gk{x-y)(t){y)ds, 
Jdn2 
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where Gk is Faddeev's Green function defined by 

If e*^"*^ 



Once equation (p!9|) had been solved, we would substitute the result into 

(20) t(k)= [ e'~^^{A^-Ai)il;{-,k)ds, 

where t is called the scattering transform, and Ai is the DN map for the unit 
conductivity. For each fixed x G 17, we would solve the following integral formulation 
of the D-bar equation: 

(21) /i(x, A:) = 1 + -1- / ^i^e^(^'-+^^)M^^^l^^^; 

(27r)^ [k - k')k' 

then the conductivity would be perfectly reconstructed as 7(0;) = /i(x,0)^. 

However, since our starting point in practice is the matrix L^, we need to reg- 
ularize the above ideal approach as explained in [33]. The matrices and Li we 
already have, and a matrix representation S/c for the single-layer operator Sk can be 
computed numerically, provided we have numerical evaluation routines for gk{x)^ 
see [22]. We expand e*^^|aj^2 ^ vector g in our finite trigonometric basis (|9j) and 
set 

(22) ^fc := [/ + Sfe(L^-Li)]-V. 

for k ranging in a fine grid inside the disc \k\ < i?, where the truncation radius 
R> is ideally chosen according to the size of noise. The choice of R falls outside 
the scope of this paper, so we will compute below reconstructions with R ranging 
in an interval. We define the truncated scattering transform by 



(23) tR{k) 




- Li)^k){x)ds for 1^1 < R, 
otherwise. 



where J^~^ denotes transforming from the Fourier series domain to the function 
domain. Finally we use the numerical algorithm in [30 to solve equation ([2T]) with 
t replaced by ti^ and denote the solution by iir{x^ k). Then 7(0;) ~ jiiR{x^O)'^. 

7. Numerical results 

We define several conductivity distributions a G L'^{Qi) on the unit disc fti = 
^(0,ri) = 1^(0,1) and compare reconstructions computed with and without the 
boundary correction procedure. 

Before proceeding with the examples, though, we need to choose an optimal 
radius r2 to be used in the boundary correction step. We do this by examining 
numerically the simple case of the unit conductivity a = 1 and using the procedure 
(a'-c^) explained in the introduction. The numerical parameters used in this pro- 
cedure are the same as in the example reconstructions, and they are given later in 
this chapter. 

We take N = 16 and simulate non- noisy ND map Ri using the finite element 
method with 1048576 triangles in Qi as explained in Section [2l Using the standard 
square norm for matrices, this yields efem = ||^i^ — i^i ||/||i^i^|| ~ 0.0000173, where 
Rf^ is the analytically calculated ND matrix for the unit conductivity. Furthermore, 
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we construct noisy ND map Rl with formula (fT2]) and c = 0.00001, giving 

\\Rl-Ri\\/\\Ri\\ ^ 0.0001 >efem- 

To avoid notational clashes, we denote by 

L^=i the DN map on 91^2 computed from noisy ND map using formula ([18]), 
L^=i the DN map computed from non- noisy ND map using formula (p!8]) . 
L^^]^ the DN map computed directly on 9(^2, 

where by 7 = 1 we mean the conductivity a = 1 extended by ([16]). The left plot 
in Figured] shows the behaviour of the error — L^^^ II /II ^7=1 II a function 

of r2. The condition number of the matrix L^^i — L^^ ranges between 1 and 20. 
It seems that we should choose r2 > 1.2. Further, the right plot in Figure [1] shows 
the behaviour of the error ||I/^=i — 1/7=1 1|/|| 1/7=1 1| as function of r2. The error 
decreases as r2 grows; it shows how the data measured on dCti contributes less and 
less to L^^i as r2 gets larger. This observation is in agreement with the known fact 
that in bit it is more difficult to obtain information from the deeper parts of the 
object [15]. 

Based on the above numerical investigation we choose r2 = 1.2 for the rest of 
this paper. We work with the following four example conductivities: 

• Example one: conductivity has a high contrast bump right on the bound- 
ary dCti and a circular inclusion near the boundary. All deviations from 
background conductivity 1 satisfy (j{x) > 1. 

• Example two: similar to Example one but with a larger inclusion having 
higher conductivity. 

• Example three: conductivity has high-contrast behaviour near dfli^ but the 
maximum of the deviation from background is not right at the boundary. 

• Example four: crude model of a cross-section of an industrial pipeline, 
similar to the case in [46 . There is a sediment layer on the bottom of the 
tube, and two round low-conductivity inclusions. 

See Figure [2] for plots of the example conductivities and their traces on d^i. 

We simulate noisy eit data for each example conductivity using c = 0.00001. 
The error \\R'^ - R^\\/\\R^\\ ranges between 0.00011 and 0.00076. 

We use the method of Section [3] with M = 32, = 6, a = 4 and 100 different 
angles to compute approximately reconstructed traces g on d^i. See the right 
column of Figure [2] for the result. Then, we compute the extended conductivity 
7 in the disc Q2 = ^(0,r2) = 1^(0,1.2) using ([T6]) and ([3]) with the radius rg = 
ri + 7/8(r2 — ri) = 1.175. Since g is only approximately the same as a\QQ^ there 
are discontinuities in 7 in all cases. 

We compute the intermediate DN maps A*-^ using the finite element method and 
425984 triangles in the annulus \ ^i- To check the accuracy of formula (flS)) 
we also calculate A^ directly by the finite element method (and 1081344 triangles 
in (^2) using the knowledge of 7. The error \\L^ — where is the DN 

map calculated directly on the boundary 9(^2, was less than 2.2% in all cases. The 
condition number of the matrix L^. — L^^ used in ([T8]) was less than 27 in all test 
cases. 

Figure [3] illustrates how the noise and the boundary correction procedure affect 
the scattering transform in example four. The first row shows the the real and 
imaginary parts of ([23]) substituting in place of L^. The second row shows the 
same functions using LJ, and the third row is again the same, but uses calculated 
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from (p!8|) . The real part of tR{k) is in the left column, the imaginary part on the 
right. The scattering transform is calculated in a grid of spectral parameters /c, 
where \k\ < 10. In white areas we have |ti^(/c)| > 15, meaning the calculation has 
failed or is close to failing due to computational error caused by large values of \k\. 

For all truncation radii R = 3.0, 3.2, 5.8, 6.0, as explained in Section [6l a 
reconstruction is calculated with and without the boundary correction procedure 
using the same reconstruction points. The conductivities and their extensions are 
pictured in figure [2l Full error graph showing -error for every reconstruction is 
pictured in figure SI Reconstructions and the corresponding errors are pictured in 
figures [SlSTTl and [51 The first reconstruction pair is always calculated with = 3, the 
second one is the one with the lowest numerical -error for the boundary corrected 
reconstruction, and the third one is with R = 6 to show how the reconstructions 
fail. 

8. Conclusion 

Our aim in this work is to find examples of simulated conductivities that (i) share 
features of conductivities appearing in applications of electrical impedance tomog- 
raphy, and (ii) allow higher-quality reconstruction when boundary correction step 
is added to the D-bar method. After experimenting with a large number of candi- 
date conductivities we concluded that for conductivities which behave moderately 
at and near 91^1, the method of approximating the trace of conductivity by an 
optimal constant is good enough. More precisely, the errors caused by measure- 
ment noise in Steps (a^) and (c^) prevented the boundary correction procedure from 
enhancing the reconstructions. 

However, we were able to find several examples where the boundary corrected 
D-bar method does provide better imaging quality than the non-corrected method 
both in terms of quantitative error and visual inspection. Four such examples are 
presented in Section [71 and all of them have high contrast features in the conduc- 
tivity right at the boundary. Consequently, most medical applications do not need 
the boundary correction procedure, but it may be beneficial or even necessary for 
some nondestructive testing, industrial process monitoring or geophysical sensing 
applications. 
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Figure 1. Left: relative error \\Ll^^i - L'^^-^\\/\\L'^^-^\\ as a func- 
tion of r2. Here i^^^i is the DN map on 80.2 computed from noisy 
ND map using formula (fT8|) and is the DN map calculated 

directly on dVL2. Here || • || denotes the standard square norm for 
matrices. By 7 = 1 we mean a = 1 extended by (fT6|) . Right: rel- 
ative error — i^7=i ||/||i^7=i|| as a function of r2. Here L^=i 
is the DN map computed from non-noisy ND map using formula 
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Figure 2. Left column: Example conductivities 7 shown in the 
extended domain (^2, the white circle indicates the inner boundary 
dVti. Right column: actual traces of the conductivities at the inner 
boundary dVti (solid line), and approximate traces at dVti (dashed 
line) whose reconstruction is explained in Section [3l 
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Figure 3. The scattering transform in example four. The first 
row shows the the real and imaginary parts of (|23|) substituting L^r 
in place of Ci% . The second row shows the same functions using L J , 
and the third row is again the same, but uses calculated from 
([18]). The real part of tji{k) is in the left column, the imaginary 
part on the right. The scattering transform is calculated in a grid 
of spectral parameters /c, where \k\ < 10. In white areas we have 
|ti?(/c) I > 15, meaning the calculation has failed or is close to failing 
due to computational error caused by large values of \k\. 
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Figure 4. L^-error graphs as a function of truncation radius R 
of the scattering transform for different examples; sohd hne is for 
the traditional D-bar reconstructions, dashed line is for boundary 
corrected reconstructions. The R -axis is the same in all four plots. 
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Figure 5. Example 1 reconstructions; the original conductivity 
on the left, traditional D-bar reconstructions on the upper row and 
boundary corrected reconstructions on the lower row; the numbers 
beside the pictures are -errors, for the full error graph, see figure 
m The first reconstruction pair is always calculated with = 3, 
the second one is the one with the lowest numerical -error for 
the boundary corrected reconstruction, and the third one is with 
= 6 to show how the reconstructions fail. 
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Figure 6. Example 2 reconstructions; the original conductivity 
on the left, traditional D-bar reconstructions on the upper row and 
boundary corrected reconstructions on the lower row; the numbers 
beside the pictures are -errors, for the full error graph, see figure 
m The first reconstruction pair is always calculated with = 3, 
the second one is the one with the lowest numerical -error for 
the boundary corrected reconstruction, and the third one is with 
= 6 to show how the reconstructions fail. 
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Figure 7. Example 3 reconstructions; the original conductivity 
on the left, traditional D-bar reconstructions on the upper row and 
boundary corrected reconstructions on the lower row; the numbers 
beside the pictures are -errors, for the full error graph, see figure 
m The first reconstruction pair is always calculated with = 3, 
the second one is the one with the lowest numerical jC^ -error for 
the boundary corrected reconstruction, and the third one is with 
= 6 to show how the reconstructions fail. 
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Figure 8. Example 4 reconstructions; the original conductivity 
on the left, traditional D-bar reconstructions on the upper row and 
boundary corrected reconstructions on the lower row; the numbers 
beside the pictures are -errors, for the full error graph, see figure 
m The first reconstruction pair is always calculated with = 3, 
the second one is the one with the lowest numerical jC^ -error for 
the boundary corrected reconstruction, and the third one is with 
= 6 to show how the reconstructions fail. 



